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Abstract 

-  A  computer  analysis  with  user  guidelines  to  analyze 
partially  continuous  multi-span  beams  is  presented. 
Partial  continuity  is  due  to  rotational  slip  which  occurs 
at  spliced  joints  at  the  supports  of  continuous  beams 
such  as  floor  joists.  Beam  properties,  loads,  and  joint 
slip  are  input;  internal  forces,  reactions,  and  deflec¬ 
tions  are  output.. 


Introduction 

Floor  joist  deflection  and  maximum  design  moment  are 
decreased  by  utilizing  continuous  joists  over  two  or 
more  spans  compared  to  simply  supported  joists. 
Splices  are  used  to  attain  the  lengths  required  for  con¬ 
tinuous  members.  The  splices,  whether  nailed,  bolted, 
glued,  or  truss  plated,  if  located  at  an  interior  support, 
will  have  rotational  slip  occur  when  subjected  to  bend¬ 
ing  moment.  Thus  the  beam  continuity  is  disrupted  with 
the  joist  acting  somewhere  between  a  simply  supported 
and  a  fully  continuous  beam. 


allows  the  input  of  some  percentage  of  complete  fixity 
but  has  the  disadvantage  of  not  being  able  to  input  the 
specific  amount  of  joint  slip  since  its  relationship  to 
the  fictitious  member  stiffness  is  not  known. 

This  study  presents  a  computer  analysis  method  and 
input  user  guidelines  to  determine  internal  forces,  reac¬ 
tions  and  deflections  of  continuous  beams  with  rota¬ 
tional  slip  at  supports.  Although  developed  specifically 
for  two-span  floor  joist  analysis  and  design,  the  method 
and  computer  program  are  applicable  to  any  con¬ 
tinuous  beam  structure. 

Joist  design  is  based  on  satisfying  both  stiffness  and 
strength  requirements.  The  stiffness  criterion  most 
commonly  used  Is  to  limit  maximum  joist  deflection  in 
a  floor  system  to  span/360  when  subjected  to  static  live 
load  of  40  pounds  per  square  foot.  The  strength 
criterion  limits  bending  stress  to  an  allowable  value 
based  on  species  and  grade  of  lumber.  Design  tables 
( 3 f  are  available  for  floor  joist  design  for  either  simple 
or  two-span  continuous  beams.  The  continuous  beam 
tables  assume  full  moment  capacity  over  the  entire 
span  (l.e.,  no  splices  or  other  loss  of  continuity). 


The  finite  element  method  of  analysis,  using  discrete 
elements  referred  to  as  matrix  structural  analysis,  Is 
the  state-of-the-art  method  used  to  analyze  continuous 
structures.  The  usual  procedure  in  matrix  structural 
analysis  to  account  for  partial  fixity  is  to  model  a  short 
or  fictitious  member  with  a  low  stiffness  value.  This 
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The  National  Association  of  Homebuilders  (NAHB)  (2) 
and  American  Plywood  Association  (APA)  (1)  in¬ 
vestigated  two-span  continuous  beams  with  splices 
near  the  inflection  points.  The  splices  were  designed  to 
transmit  shear  force  since  moments  near  the  inflection 
points  are  assumed  small.  Splice  slip  was  not  con¬ 
sidered. 


>  Italicized  numbers  In  parenthesis  refer  to  literature  cited  at  the  end  of 
this  report. 


Splice  slip  Is  due  to  the  behavior  of  the  mechanical 
fastening  system.  Nails  (and  truss  plates)  slip  due  to 
bending  of  the  nail  (tooth)  and  crushing  of  the  wood. 
Many  studies  have  established  nonlinear  load-slip 
behavior  for  nailed  joints  In  single  shear.  When  bolts 
are  used  they  are  installed  in  oversized  holes;  thus  a 
certain  amount  of  movement  occurs  before  loads  are 
transmitted.  Slip  measurements  related  to  design  load 
levels  are  required  for  various  splice  configurations. 

Splice  slip  affects  both  joist  stiffness  and  strength. 
Slippage  of  a  splice  at  an  interior  support  will  result  in 
increased  deflection  throughout  the  span  and  in  a 
decreased  moment  capacity  at  the  support. 

Matrix  analysis  (4)  of  continuous  beams  assumes  sup¬ 
ports  either  fixed  or  pinned.  Intermediate  end  condi¬ 
tions  are  generally  modeled  by  addition  of  a  fictitious 
member  with  low  stiffness  or  a  rotational  spring.  Matrix 
analysis  assumes  each  joint  having  discrete 
displacements  corresponding  to  degrees  of  freedor  ; 
thus  no  discontinuity  in  displacement  can  occur.  Loads 
at  or  between  joints  and  support  settlements  can  be  in¬ 
cluded  in  the  analysis. 

Theory 

The  philosophy  of  the  analysis  is  illustrated  in  figure  1. 
Deflections  due  to  loads  are  found  by  matrix  method 
assuming  no  joint  slip.  An  experimentally  determined 
slip,  0S,  is  input  and  allocated,  dr.  0j,  to  the  adjacent 
right  and  and  left  spans.  Deflections  caused  by  adja¬ 
cent  slips  are  calculated  and  superimposed  on  the 
deflections  due  to  loads.  Member  end  actions  and  sup¬ 
port  reactions  are  calculated  from  the  member  deflec¬ 
tions. 

The  member  stiffness  matrix  assumes  the  member  be¬ 
ing  subjected  to  lateral  loads  which  induce  bending 
moments  and  shear  forces.  The  effects  of  shear  are 
neglected;  thus  two  degrees  of  freedom  at  each  node, 


vertical  and  rotational,  are  necessary.  The  degrees  of 
freedom  are  numbered  sequentially  from  left  to  right; 
thus  the  vertical  translation  of  the  left  end  is  1,  the 
rotation  2,  and  those  of  the  right  end  are  3  and  4, 
respectively,  as  shown  in  figure  2.  Also  shown  is  the 
member  stiffness  matrix  with  forces  corresponding  to 
unit  displacements  for  each  degree  of  freedom.  The 
positive  sign  convention  for  the  member  forces  and 
displacements  are  as  shown,  with  vertical  translation 
upward  and  rotation  counterclockwise  being  positive. 

The  beam  is  modeled  with  nodes  at  points  of  support, 
changes  in  cross-section  properties,  and  at 'any  other 
point(s)  where  shear  forces  and  bending  moments, 
and/or  vertical  and  rotational  displacements  are  to  be 
computed.  Figure  3  shows  node  locations  for  an  exam¬ 
ple  two-span  beam.  It  includes  the  support  locations 
plus  an  arbitrary  interior  location,  a  distance  x  from  the 
left  support,  where  forces  and  displacements  are  to  be 
calculated.  Again,  there  are  two  degrees  of  freedom  at 
each  node;  they  are  numbered  sequentially  starting  at 
the  left  end  of  the  beam.  Member  numbers  (circled)  are 
also  sequential  from  the  left  end.  Positive  node  forces 
and  displacements  are,  as  shown,  upward  and 
counterclockwise.  The  member  stiffness  matrices  are 
superimposed,  as  shown  symbolically,  to  form  the 
structure  stiffness  matrix,  [Sj. 


The  structure  stiffness  matrix  is  rearranged  and  parti¬ 
tioned  related  to  the  unknown  displacements,  Dg,  and 
the  known  boundary  displacements,  Ds,  representing 
support  conditions: 


[S]  = 


sdd  i  Sds 
Ssd  1  Sss 


where  the  subscript  s  refers  to  support  degrees  of 
freedom  (with  known  displacements  of  zero  or  support 
settlement  value)  and  d  refers  to  degrees  of  freedom 
with  unknown  displacements. 
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Figure  1 Superposition  of  load  and  slip  deflected  shapes. 
(M  149  237) 


2 


2(  im  member 

f7  V 

K -  L  - 


From  figure  4b,  the  moment  equilibrium  equation  at  the 
interior  support  is: 


Mj  =  -  Mr 


(2) 


(a)  Member 

F I2_EI  6EI 

-I2EI 

6EI 

where  Mj,  Mr  are  the  moments  corresponding  to  Of  and 
0r.  From  standard  beam  theory,  the  relationship  be¬ 
tween  moment  and  rotation  is: 
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where  E  is  the  modulus  of  elasticity,  1  the  moment  of 
inertia,  and  L  the  length  of  the  left  and  right  spans 

6EI 

2EI 

-6EI 

4EI 

(subscripts  i,  r),  respectively.  The  relationship  is  ap¬ 

L2 

L 

L2 

L 

plicable  to  members  having  a  rotation  at  one  end  of  the 
member. 

(b)  Member  Stiffness  Matrix 
Figure  2  — Member  degrees  ot  freedom  and  stiffness  matrix. 
(M  149  238) 


Node  forces  and/or  support  settlements  (if  any)  are  in¬ 
put  with  sign  convention  in  figure  3b.  Fixed  end  shears 
and  moments  corresponding  to  member  loads  are  input 
with  sign  convention  in  figure  2a. 
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(a)  real  beam 


The  unknown  joint  displacements,  {D<j},  due  to  load, 
{a},  are  found  by: 

M*M  "W 

IA(  =  jM  -  |Afemj  ~  [Sds]  )■>»{ 

where  {A}  is  the  general  load  matrix  consisting  of  the 
specified  joint  loads,  {Aj\,  fixed  end  reactions  due  to 
between-the-joint  loads,  (Afem},  and  the  forces, 
[Sdsj{Ds},  due  to  support  settlement,  {Ds}. 

The  known  joint  slip  is  allocated  to  adjacent  spans 
based  on  compatibility  and  equilibrium  as  indicated  in 
figure  4  in  which  moments  and  rotations  are  shown  in 
the  positive  direction.  The  effects  of  spans  other  than 
the  adjacent  spans  are  neglected.  The  sign  convention 
for  the  slip  rotation  is  that  the  slip  angle  is  measured 
from  the  tangent  to  the  elastic  curve  in  the  right  span 
to  the  tangent  of  the  elastic  curve  in  the  left  span  with 
counterclockwise  rotation  being  positive. 

From  figure  4a,  the  compatibility  equation  is: 

0S  =  0,  -  8r  (1) 

where  6 s  is  the  slip  rotation,  and  Of,  0r  are  the  rota¬ 
tional  allocation  of  the  slip  to  the  left  and  right  adja¬ 
cent  spans,  respectively. 


(b)  MODELED  BEAM 


(C)  STRUCTURE  STIFFNESS  MATRIX 


Figure  3.— Structure  degrees  ot  freedom  end 

stiffness  matrix.  (M  149  239) 
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tb)  EQUILIBRIUM 

Figure  4.— Joint  slip  compatibility-equilibrium  relations. 
(M  149  240) 


Substituting  equations  (3)  into  (2)  yields  a  relationship 
between  and  0r,  and  substituting  this  into  equation 
(1)  results  in  the  following  allocation  of  slip  to  adjacent 
spans  (assuming  Er  =  E|): 


(4a) 


(4b) 

To  determine  displacements  due  to  a  known  slip,  the 
continuous  beam  is  separated  into  two  structures  at 
the  support  where  the  slip  occurs.  The  slip  allocations, 
9f  and  0r,  are  applied  as  support  displacements  to  each 
side  of  the  separated  structure.  Deformations  are  then 
determined  by  matrix  analysis. 

The  member  end  displacements  due  to  slip  are 
superimposed  on  the  displacements  due  to  loads. 
Member  end  displacements  multiplied  by  the  member 
stiffness  matrix  summed  with  the  fixed  end  moments 
due  to  load  result  in  member  end  forces.  Adjacent 
member  end  forces  combined  with  applied  joint  forces 
determine  support  reactions. 

Procedure 

The  Fortran  program  using  this  theory  is  given  in  ap¬ 
pendix  A.  The  program  is  intended  to  be  as  complete 
as  possible  so  that  It  can  be  easily  modified  for  future 
research;  It  is  not  a  production  tool  and  no  effort  has 
been  made  to  make  it  as  efficient  as  possible.  The 
units  used  must  be  compatible;  the  example  in  this 
paper  uses  Input  lengths  In  Inches,  forces  In  pounds, 
moments  In  inch-pounds,  slip  in  radians,  modulus  of 
elasticity  in  pounds  per  square  Inch,  and  moment  of  in¬ 
ertia  in  Inches'.  Output  is  in  the  same  units. 
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Figure  5.— Schematic  representation  ot  the  makeup  of 
computer  deck.  (M  149  241) 


The  program  is  limited  to  one  slip  per  member,  and  no 
slip  at  either  the  first  or  last  support.  Slips  at  adjacent 
joints  are  permissible  provided  each  is  indexed  to  a 
separate  adjacent  member.  The  program  Is  arbitrarily 
limited  to  a  maximum  of  10  elements  which  is  con¬ 
sidered  adequate  for  most  2-  or  3-span  floor  joist 
systems  (extra  nodes  may  be  included  between  sup¬ 
ports).  For  larger  problems,  this  limitation  can  be 
removed  by  modification  of  the  dimension  statement. 

The  degrees  of  freedom  are  numbered  sequentially 
from  the  leftmost  support  with  the  vertical  translation 
first  and  the  rotation  second.  Sign  conventions  have 
been  previously  defined  and  are  shown  in  the  positive 
direction  in  figures  2,  3,  and  4. 

The  Fortran  program  is  stored  on  both  tape  and 
punched  cards  on  the  FPL/MACC  system.  Input  data  re¬ 
quired  are  described  in  table  1.  The  command  se¬ 
quence  for  the  FPL/MACC  system  to  access  the  tape  is 
given  In  table  2;  that  for  the  punched  deck  is  shown  in 
figure  S. 

The  following  example  illustrates  both  input  data  re¬ 
quired  and  output  generated. 

Example 

The  two-span  beam  of  figure  3a  is  assumed  to  have  a 
splice  at  B  which  slips  +  0.00432  radian  when  loaded. 
Other  values  for  this  example  are: 

W,  =  SOIblft  =  4.16  lb/ln.  I,  -  I,  -  20.8  In* 

L,  *  12  ft  -  0  in.  *  144  In.  W,  a  10  lb/ft  *  0.833  Ibfln. 

E,  a  E,  *  1,700,000  lb/ln.'  1.,  a  9  ft  -  0  In.  a  108  In. 

Shear  and  moment  diagrams,  and  deformed  shape  for 
this  partially  continuous  beam  example  are  required. 
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Tabic  1.— Data  card* 


Number  ol 
cards 
required’ 

Information  required1 

Input  data 

In  columns 
numbered 

Fortran 

format 

1 

a.  NM  =  Total  number  of  members 

1  through  3* 

13 

b.  NS  =  Total  number  of  supports 

4  through  6* 

13 

corresponding  to  degree  of 
freedom  numbers  (i.e., 

NS  a  2  for  shear  and 
moment  at  fixed  support) 
c.  NF  =  Total  number  of  joint  slips 

7  through  9* 

13 

d.  NA  =  Total  number  of  degrees  of 

10  through  12‘ 

13 

NS 

freedom  corresponding  to 
joint  loads  (l.e.,  NA  =  2 
for  joint  with  applied 
vertical  load  and  moment) 

For  each  card  (support): 

a.  Structure  degree  of  freedom  number 

1  through  3‘ 

13 

b.  Support  settlement  (in.  or  radians) 

4  through  13 

F10.6 

NM 

For  each  card  (member): 
a.  Member  number 

1  through  3* 

13 

b.  Modulus  of  elasticity  (lb/in.3 

4  through  13 

F10.0 

c.  Moment  of  inertia  (in.4) 

14  through  21 

F8.2 

d.  Length  (In.) 

22  through  28 

F7.2 

e.  Left  end  fixed  end  shear  due  to 

29  through  37 

F9.2 

member  toads  (lb) 
f.  Left  end  fixed  end  moment  (in.-lb) 

38  through  46 

F9.2 

g.  Right  end  fixed  end  shear  (lb) 

47  through  55 

F9.2 

h.  Right  end  fixed  end  moment  (in.-lb) 

56  through  64 

F9.2 

NA 

For  each  card  (joint  load): 
a.  Structure  degree  of  freedom 

1  through  3* 

13 

number  corresponding  to 
joint  load 

b.  Joint  load  (lb  or  in.-lb) 

4  through  13 

F10.2 

NF 

a.  Member  number  (may  be  either 

1  through  3* 

13 

left  or  right  span  adjacent  to 
slip) 

b.  Member  degree  of  freedom  number 

4  through  6* 

13 

(either  2  or  4  corresponding  to 
member  selected  in  “a") 
c.  Structure  degree  of  freedom 

7  through  9* 

13 

number 
d.  Slip  (radians) 

10  through  19 

F10.6 

'  Cards  must  be  sequenced  in  this  order. 
1  Zero  values  may  be  entered  by  blanks. 

*  Values  must  be  right-adjusted. 
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Table  2.— Tape  command  sequence 

©RUN... 

©PASS... 

©CAT  HUD*CONTINUBEAM. 

©ASG,  AX  HUD'CONTINUBEAM. 

©ASG,  TH  DIMEN  ‘LUMBER., U9H,  7639 
©TGET  DIMEN ‘LUMBER.,  HUD'CONTINUBEAM. 
©XQT  HUD'CONTINUBEAM.SLIP 
Data  as  per  table  1 
©FIN 


The  beam  is  modeled  with  joints  at  points  of  support 
and  at  other  arbitrary  locations  deemed  necessary  to 
define  the  shear  and  moment  diagrams  and  the  de¬ 
formed  shape.  For  illustrative  purposes,  only  one  ar¬ 
bitrary  location  5  feet  (60  in.)  from  the  left  support  (x  = 
60  in.)  is  selected.  Thus  the  beam  is  modeled  with 
structure  degrees  of  freedom  and  members  numbered 
sequentially  from  the  left  as  in  figure  3b. 

The  input  data  corresponding  to  table  1  are: 

Card  1 


NM  =  3 

NS  =  3  (corresponds  to  structure  degrees  of 
freedom  1,  5,  and  7) 

NF  =  1  (corresponds  to  structure  degree  of  freedom 
6) 

NA  -  0  (no  joint  loads,  only  member  loads) 

Cards  2,  3,  and  4,  respectively 
Structure  degree  of  freedom  =  1,5,  and  7, 
respectively 

Support  settlement  =  0  (all  cards) 


Cards  5,  6,  and  7,  respectively 
Member  number  =  1,2,  and  3,  respectively 
Modulus  of  elasticity  =  1,700,000  lb/in.]  (all  cards) 
Moment  of  inertia  =  20.6  in.4  (all  cards) 

Length  =  60,  84,  and  108  in.,  respectively 
Left  end  fixed  end  shear  for  member  1  = 


W1X  50x5 
2  2 


+  125.01b 


Left  end  fixed  end  moment  for  member  1  = 


W1X’  50  x  5J  x  12 

_ 


+  1250.0  in.-lb 


Fixed  end  reactions  for  right  end  and  for  other 
members  similarly  found. 

(Note,  sign  convention  as  per  figure  2;  thus  right  end 
fixed  end  moment  is  negative.) 


displacements  to  the  left  and  right  of  the  joint  are 
given. 

The  results  are  plotted  as  case  II  in  figure  6.  Results  for 
simply  supported  and  fully  continuous  beams,  found  by 
standard  structural  analysis,  are  given  as  cases  I  and 
III  for  comparison.  As  expected,  the  behavior  of  the  par¬ 
tially  continuous  beam  is  bounded  by  the  simple  and 
fully  continuous  cases.  The  displacement  at  x  =  5  feet 
is  0.6318  and  0.4126  inch  for  the  simply  supported,  and 
fully  continuous  beams,  respectively.  The  partially 
continuous  beam  results  in  a  deflection  of  0.4738  inch. 
Use  of  a  joint  at  B  reduces  the  simply  supported  deflec¬ 
tion  by  about  25  percent.  The  negative  moment  over  the 
center  support  is  reduced  from  the  fully  continuous 
value  of  558  to  406  foot-pounds  for  the  partially  con¬ 
tinuous  beam;  however, the  corresponding  positive  mo¬ 
ment  at  x  =  5  feet  is  increased  from  643  to  706  foot¬ 
pounds.  This  is  still  about  20  percent  less  than  the 
simply  supported  positive  moment.  Figure  7  details  the 
discontinuity  of  the  deformed  shape  at  support  B. 

Program  Alteration 

The  program  is  arbitrarily  limited  to  10  elements  (NM  = 
10);  the  corresponding  number  of  joint  degrees  of 
freedom  is  22  (2NM  +  2).  Dimensioned  arrays  have 
values  of  10,  22,  or  4  corresponding  to  number  of 
elements,  number  of  joint  degrees  of  freedom  or 
number  of  degrees  of  freedom  per  element.  To  increase 
the  number  of  elements,  change  the  dimension 
statements  in  the  first  six  cards  of  the  source  program 
as  follows: 

a.  Dimensions  of  value  10  are  increased  to  the  new 
number  of  elements,  NM. 

b.  Dimensions  of  value  22  are  increased  to  a  value  = 
2NM  +  2. 

c.  Dimensions  of  value  4  are  unchanged. 

d.  Dimensions  of  the  array  called  SCRACH  are  in¬ 
creased  to  3  x  (2NM  +  2). 


No  NA  cards  (NA  =  0) 

Card  8  (NF  =  1) 

Member  number  of  left  adjacent  span  =  2  (alternate¬ 
ly,  member  3,  the  right  adjacent  span  could 
be  selected) 

Member  degree  of  freedom  (figure  2a)  =  4  (alternate¬ 
ly,  member  degree  of  freedom  =  2  could  be  selected 
to  correspond  to  member  3)  Structure  degree  of 
freedom  (figure  3b)  =  6  Slip  =  +  0.00432  radian 

The  output  forces  and  displacements  are  given  in  ap¬ 
pendix  B  to  illustrate  the  output  format  (which  is 
referenced  to  the  structure  degree  of  freedom  number¬ 
ing).  Joint  displacements  are  not  given  directly  since 
the  slip  creates  a  discontinuity;  member  end 
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-  CASE  I  -- 

-  CASE  n  -- 

-  CASE  m  * 


SIMPLY  SUPPORTED  BEAMS 
PARTIALLY  CONTINUOUS  BEAM 
CONTINUOUS  BEAM 


W-  50  lbs/ ft 


W  -  !Q  Ibs/ft. 


a 


t 


!2  ’-o" _ 

E  --  /,  700,000  psi 
I*  21  in? 


g'.  o" 

E  -- 1,700,000  psi 
1=2!  in? 


Rj  -  300  lbs 
Ra-  266.2  lbs 
Ra=  253.5  lbs 


Rj  --  345  lbs 
Ra*  423.9  tbs 
Rm -  453.4  lbs 


Rj  *  45  lbs 
Ra=  -0.1  lb 
Rm -*  -17.0  lbs 


SHEAR 

(tbs) 


MOMENT 
(ft.  lbs) 


DEFORMED 

SHAPE 

(in.) 


CASE  I  CASE  a  case  m 


Figurs  7. — Exsmplt  tpllc*  slip  dstsll  at  support. 


(U  149  316) 


r 


Appendix  A 

Computer  Source  Program 

1  C  MATRIX  ANALYSIS  OF  CONTINUOUS  BEAMS 

2  C  MAXIMUM  NUMBER  OF  SPANS  =  10 

3  DIMENSION  JNS (22) . DSC22) , JNSSC22) . INSC22) . AJ C22) .SRC22) .DDDC22) ,  T 

4  1C22),BC10),V(10),UC1B.),XC10),SM(10,4.4},EC10),XIC10),XLC10),AM(10, 

5  24) ,  AMF  ( 10,4),COORC22),XLSC  101  .SC22.22)  .A  (22).  JNDC22)  .RRC22.22)  ,R(2 

6  32.22) .SDDC20.20) , SSDC22.22) .SDSC22.22) .SSSC22.22) , ASC22) ,ADC22) ,CC 

?  422),YC22).SDDI (20.20). SCRACH (60.60). DD (22), D (22), SLIP (22). SLIP  1(22 

8  5 ) ,  SL  IP2  C 22 ) ,  DL  ( 10 . 4) .  DR  C 1 0 , 4) , AMM C 10 , 4) ,  A  J  J  C  22 ) 

9  C 

10  C  INPUT  TATA 

11  READ  (5,  DNM.NS.NF.NA 

12  DO  5  1=1, NS 

13  READ (5, 2) JNS ( I) .DSC  I) 

14  5  JNSSC  I)  -JNSCI) 

15  1  FORMAT (413) 

16  2  FORMAT ( 13, F 10 . 6) 

IT  C 

18  C  INITIALIZE  CONDITIONS 

19  N IF=0 . 

20  NSSJ J=0 . 

2 1  NMM=NM 

22  NSS=NS 

23  DO  10  1=1,22 

24  INS (I >  =0. 

25  AJJ ( I ) =0 . 

26  AJ ( I) =0 . 

2Z  SR(I)=0. 

28  10  DDD ( I ) =0 . 

29  C 

30  C  MEMBER  STIFFNESS  MATRIX 

31  DC  15  J  =  1 , NM 

32  READ  (5. 3)  M,  E  (M) ,  XI  (M) ,  XL  (M) ,  (AM(M,  I),  1  =  1,4) 

33  3  FORMAT C 13, F 10 . 0, F3 . 2.F7 . 2, 4F9 . 2) 

34  B  C  M)  =6 .  *E  C M)  *X I  CM)  /XL  CM)  **2 . 

35  V(M) =2 . *B  CM) /XL (M) 

36  UCM) =4.*ECM)*XICM)/XLCM) 

3?  X(M) =U(M)/2. 

38  SM(M, 1 . 1 ) =V(M) 

39  SM(M, 2, 2) =U(M) 

40  SMCM.3, 3) =V(M) 

41  SM(M,  4. 4)  =(J(M) 

42  SM(M, 1,3) =-V(M) 

43  SMCM.3, 1)=-V(M) 

44  SMCM.2,4) =X(M) 

45  SMCM, 4,2) =X(M) 

46  SMCM, 1,2) =B(M) 

4?  SMCM, 1,4) =B CM) 

48  SMCM,  2, 1 ) =B  CM) 

49  SMCM,  4,  1 )  =B  CM) 

50  SMCM, 3, 2) *-B  CM) 

51  SMCM. 3,4) =-B CM) 

52  SMCM, 2, 3) *-BCM) 

53  15  SMCM. 4,3) =-BCM) 

54  DO  16  J  = 1 ,NMM 

55  DO  16  L=1 .4 

56  16  AMF C J, L) *AMC J , L) 

5?  C 

58  C  STRUCTURE  SPAN  LENGTHS 

59  NJ*2*NM+2 
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btl 

NO J =N  J 

b  1 

co  or  f  n=o. 

62 

DO  20  I=3,NJ,2 

63 

(  =  u  -  n -2 

64 

20 

COOR ( I ) =C00R ( 1-2) +XL (K) 

65 

NSL =NS- 1 

66 

1  =  1 

6 17 

DO  35  J=3,UJ,2 

68 

no  30  K=2,.  NS 

69 

IFU.EQ.  JNSCK)  )G0  TO  25 

70 

GO  TO  30 

71 

25 

JJ=JNS (K- 1 ) 

7  2 

XLS C  I ) =C00R ( J) -COOR  C JJ) 

73 

1  =  1  +  1 

!  -4 

GO  TO  35 

■7tr 

30 

CONTINUE 

76 

77  C 

7C 

b-J 

CONTI HUE 

78  C 

STRUCTURE  JOINT  LOADS 

79 

IF (NA) 50. 50, 40 

80 

40 

DO  45  1  =  1..  NA 

81 

READ  <  5. 4)  NAM..  AJ  (NAM) 

82 

4 

FORMAT  (.  13, F  10.2) 

8  b 

A JJ (NAM) =AJ (NAM) 

84 

45 

CONTINUE 

85 

86  C 

50 

CONTINUE 

o  r' 

999 

CONTINUE 

88  C 

89  C  STRUCTURE  STIFFNESS  AND  LORD  MATRIX 


90 

DO  55  1=1,22 

91 

AC  I ) =0. 

92 

DO  55  J  =  1 , 22 

93 

55 

S  C I ,  J ) =0 . 

94 

DO  75  1=1, NM 

95 

M=  I 

9b 

Ml  =201-  1 

93 

M2 =2  *11 

90 

113=2+11+1 

99 

114=2+11+2 

100 

ri=n+HSSjj/2 

101 

SC  11 1,113)  =  SMC  11,  1,3) 

102 

sail,  114)  =SI1CM,  1,4) 

103 

S  CM2 , 113)  =511(11,2,3) 

104 

S  (112, 114)  =51141,2,4) 

105 

SCM3,MD=SMCM,3, 1) 

1 06 

S  CM3, 112)  =SM(M,  3, 2) 

107 

S  (114,  Ml )  =SM(M,  4,  1) 

108 

SC  114,  M2)  =SMCM,  4, 2) 

109 

S (Ml , Ml ) =S(M1 , Ml ) +SMCM,  1.  1) 

1  10 

S C 11 1 , M2 )  =S  (Ml ,  112 )  +S11  CM,  1,2) 

1  1  1 

S  C 11? .  11 1 )  =S  C 112 , 11 1 )  +SMC M,  2,1) 

1 12 

S  C 112 ,  M2 )  =S  C  M2 , 112 )  +SMC  M,  2 , 2  ) 

1  13 

IF  (I -Nil)  60, 65 .70 

1  14 

60 

S  <  113 , 113 )  =S  C 113 ,  M3 )  +SMC  M,  3 , 3) 

1  15 

S  C 113 . 114 )  =  S  C 113 , 114 )  +S M  C  M,  3 , 4) 

1  16 

S  C 114 , 113 )  =  S  ( M4 , 113 )  +S11 C 11 ,4.3) 

1 17 

S  C 114, 114)  =S  (M4, 114)  +SM C 11 .4,4) 

1 18 

GO  TO  70 

1 19 

65 

S  CM3.  M3)  =  SMC  11, 3, 3) 

120 

S (113, 114)  =SM41.  3.4' 

121 

S  414,113)  =SM(11, 4. 3) 

122  SCM4.M4) =SMCM,4,4> 

123  70  CONTINUE 

124  A  ail )  =0  (  Ml )  -filial.  1 1  +flj  (Ml I 

1 25  fi(  M2 )  -fi  ( M2 )  -fill '  M.  2 )  +fi  J  C  M2 ) 

126  fi  ( M3  I  = A  ( M3 )  -filial.  3)  +fiJ  (M3I 

127  A  CM4) =fi (M4) -AM(M. 4) +AJ (M4) 

128  75  CONTINUE 

129  C 

130  C  REARRANGE  STIFFNESS  MATRIX  FOR  JOINT  RESTRAINT 

131  C  DETERMINE  JND 

132  I  -1 

133  DO  95  1  =  1. N  J 

134  DO  90  K - 1 . NS 

135  IF ( I-JNS(iO  180.95. 80 

136  80  CONTINUE 

137  I F i K-NS ) 90 . 85 . 90 

1313  85  JND  (LI  =  1 

139  L=L+1 

140  90  CONTINUE 

141  95  CONTINUE 

142  C  MOVE  '-ERTICAL  COLUMNS 

143  ND=NJ-NS 

144  DO  100  1=1. NJ 

145  DO  100  J= 1 » ND 

146  KK  =  JND  ( J) 

147  100  PRC  I , J) =SCI,KK) 

148  DO  105  1=1. NJ 

149  DO  105  J  =  1 ,  NS 

150  KK  =JNS ( J 1 

151  1 05  PR  < I . J+ND ) =S  C I . KK ) 

152  C  MOVE  HORIZONTAL  ROUS 

153  DO  110  J  =  1 . N J 

154  DO  110  1=1. ND 

155  KK =  JND ( 1 1 

156  110  Rc I, JI =RP(Kk. J) 

157  DO  115  J= 1 . NJ 

158  DO  115  1=1. NS 

159  KK= JNS (  I  > 

160  115  RCI+ND.  J)  =RR(KK.  J) 

161  C  TO  SUBDIVIDE  MATRIX 

162  DO  120  1=1. ND 

163  DO  120  J=1.ND 

164  120  SDD ( I • J) =R  C  I . J ) 

165  DO  125  1=1. NS 

166  DO  125  J= 1 . ND 

167  125  SSD 1  I . J)=P( I+ND.J) 

168  DO  130  1=1- ND 

169  DO  130  J  = 1 . NS 

170  130  SDS C I . J ) =R ( I . J+ND) 

171  DO  135  1=1, NS 

1 '2  DO  135  J= 1 , NS 

173  135  SSSC I . J) =R( I+ND. J+ND) 

174  C  REARRANGE  LOAD  MATRIX 

175  DO  140  I = 1 ,ND 

176  Y =JND 1  I  1 

17?  140  A  In  I  i  =  fi I K  i 

178  DO  145  1=1- NS 

179  t  =JNS  ( 1 1 

180  145  AS'  I  -  -  AM  i 

181  C 

182  C  TO  FIND  JOINT  DISPLACEMENTS 

183  DO  150  I  =  1.1  ID 
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184 

C  C I  ,i  =0 . 

185 

DO  150  K  = 1 , NS 

186 

150 

C  C I  i =C C I ' +SDS ( I , K) *DS ( K ) 

187 

DO  155  1=1, ND 

188 

155 

yen  =adc  n-ccn 

189 

C ALL  r  IT I NV2 ( SDD . SDD I , ND ,  20 » 20 

190 

DO  160  1=1, ND 

191 

I'D '  I)  =0. 

192 

DO  160  K  =  1 . ND 

193 

160 

I'D  f  I  '  =I'D  ( I )  +SDD  Id.  K )  *Y ( K  ' 

194 

DO  162  1=1, ND 

195 

J  =  JND  1  D+NSSJJ 

196 

DDI"  J'  =I)D(  n+DDD'J) 

197 

162 

D ( J ) =DD t  1  i 

198 

DO  165  1=1, NS 

199 

1  =  J  HE  D+NSSJJ 

20O 

DDI"!  >  =DS(  I  '+DDI'fk.) 

201 

165 

D  '  k  )  =I'S  I ) 

202  C 

203  C 

IF  SLIP  hT  SUPPORT  OCCURS 

204 

IF iNF '340,340, 170 

205 

170 

CONT I NIJE 

206 

DO  175  1  =  1.  HI  111 

207 

DO  175  J=l,4 

208 

DO  175  K =1,22 

209 

hIK  I  ,  J)  =0. 

210 

0.J  CK )  =0. 

211 

175 

DS  ( K  :>  =c. 

212 

IF (N IF-2+NF) 180.340,340 

213 

180 

N IF  =N IF+1 

214 

KN IF  =  N IF/NF 

215 

GO  TO  C  190.275) . KN IF 

216 

19P 

CONTINUE 

217 

READ (5.6) M, KKK , NFM. SL I P ( NFM) 

218 

6 

FORMAT'S  13, F 10. 6) 

219 

INS  CM) =NFM 

220 

NS  J  =2#C  (MFM+1 '  --'2)  -  1 

221 

L=0 

222 

DO  195  1=1, NSS 

d  d  i 

IFCJNSSd)  .EQ.NSJ)  GO  TO  200 

224 

I F  C  JNSS ( I ) . EQ . NFM)  GO  TO  200 

225  195  L=L  +  1 

22b  200  GO  TO (2 10. 205, 2 10, 205), KKK 

22 7  285  SLIP1  CNFM)=SLIP( N F M ) *XL S ( L ) *X  I  ( M+ 1 ) / ( X I  (1 1+ 1  > *XL S ( L )  +X I  ( M ) *XL S ( L + 1 ) 

228  1) 

229  SL IP2 (NFM) *SL IP  1 (NFM) -SL IP (NFM) 

230  210  CONTINUE 

231  C 


232  C 

FIND  NUMBER 

OF  JOINTS  AND  MEMBERS  IN 

233 

MLJ=2*C  (NFM+D/2) 

234 

HLM= (NLJ-2) 

/2 

IK 

L-  J 

236 

237  C 

NRM=NMM-NLM 

UR J =2+NRM+2 

DETERMINE  NUMBER  OF  SUPPORTS  IN  LEFT 

238 

239 

L  =0 

DO  220  1=1, 

NLJ 

240 

DO  220  J=l, 

NSS 

241 

IF  r  I - JNSS ( J 

) ) 220. 2 15, 220 

242 

215 

L  =L+ 1 

243 

220 

CONTINUE 

244 

DO  225  1=1. 

NSS 

245 

IF' JNSS'  I  » . 

EQ. NSJ) GO  TO  230 

246 

IF'  HISS'  I) . 

EO .NFM) GO  TO  230 

25  mu !  ui'JE 


14  i  230  L "L -  i 

250  -"'is  coir  ;;;:je 

?5  1  GO  1 0  .  250 , 245 . 250  •  245 ')  .  KK I 

_-.b  ‘  ""  ( 1  ;  IL5-  1  1  =2+  1  i.NFM-  1 )  ■  2  ‘  +  1 

934  2 [ : 9  .  iLr.'i  =NFM 

255  r.c  2r;,,  =SL  I P 1  (NFI1* 

550  250  C  i.’n  i !  i  HUE 

252  NT  il  ,! 

259  N2.  =  l  iLS 

2  2  0  2  U  T  Li  9  9  '9 

2b!  255  C 0 NTIH ij E 

29,2  C  DETERMINE  NUMBER  OF  SUPPORTS  IN  RIGHT  SPANS 
223  KNFM  =  2:i:i  1  NFI1-  !  '  -‘2'+l 

264  NSSj  ;=nsj-i 

22b  DO  205  I-KNFM.NOJ 

229  IF ' I-JMSSf J) >285.280.285 

229  280  !  =L  + ! 

250  253  CONTINUE 

”71  TO  290  1=1, NSS 

252  IF'  JMSST  P  .EQ.NS.J1G0  TO  295 

2  5  3  JF2  It ISS  i  I  >  .  EO  .  NFM'i  GO  TO  295 

2~1  g|  V-;  300 

220  on  TO  325 .310,325,310') ,  I  KK 
5"  5  310  NRS--L+2 

25?  J Nb  (  2 'i  =NFt1-NSSJJ 

250  JNS  •  1,'  =NFM-  1-HSSJJ 

2.2i  DS(2)  *SLIP2(NFM> 

223  DO  320  1=1. NSS 

239  I F 1 i i  iSS 1  1 1 . GT . HFM) GO  TO  315 

2  -  3  GO  TO  320 

292  3  15  LM=LI1+  1 

2-.-  JIISi  LIT  =JNSS<  D-NSSJJ 

;  320  CONTINUE 

325  CONTI MIJE 
I  if ! -f-IF'tl 
29:  I1J9P! 

392  NS  ' 1  IF 

392  GO  ; f  999 

a  ~4'"  '■ 1  ib T ;  uuE 

'  V  | ' 

TT  upep impose  deflections  and  find  member  end  forces 
MJ  545  r  1  =  1 . 1  li  t  I 

2  9-:  Ml  =  2  1 11-1 

2?b  r  12  =  I 

OfiM  M'-: -  2  1 1 1+1 


~f]  1 

2  tet¬ 

•7 

:0.' 

f'l 

ri.  i ' 

=DDD'I11 ' 

~  n  z 

DL 

i  M .  2  ' 

DDDC12  ' 

704 

DR 

>  fl.  3  > 

=DDD ' M3  J 

3Dj 

745  [.D 

I  t.  4  ' 

=DDD  '  T 14 

rob 

Du 

b-lb 

11=  1  .  NMI1 

V\7 

IT 

'  i  (F  i  7 

50.395.75i 

7  No 

350  1  - 

IIISi  r 

h 

7-  U  c* 

Ml 

-2 'll- 

1 

310 

1 12  = 

24  1 

31  1 

113  = 

2  41+1 

312 

(14  = 

2  H ! f 2 

313 

IF' 

i  .  GE  . !  1 1  .  AND  .1"  .  LE .  M4 '  GO 

TO  355 

3  14 

GO 

TO  395 

315  355 

1 1  = 

I  -241  f  2 

316 

GO 

TO  ' 395,360. 395.370) . p 

K 

3 1 7  360 

DL ' 

1 !.  2  '  =DL  '  1 1. 2  '  -SL  IF'  1  '  L) 

3  1 8 

IF  ' 

f  1-  1  '  365 . 395  -  365 

3 1 9  365 

'Tip  f. 

(1-  1 . 4)  =DP '  i  l-  1 , 4)  -SL  IP2 

(  ['  ) 

320 

GO 

TO  395 

321  370 

DP  i 

t  ( .  4  '  =DP  i  1 1, 4)  -SL  I P2  O  i 

322 

IFC 

!  I-MMM)  375 . 395  -  375 

323  325 

DL  ( 

11+1 .2  i  =DL  1 11+1  •  2  1  -  SL  IP  1 

O') 

324  395 

CONTINUE 

32  j 

DO 

400  11=1.  NMM 

326 

DO 

400  J  =  1 , 4 

327  400 

Mf  11 

If  11.  J)  =SMf  M.  J  .  1 )  *DL  (  M.  1 )  +SM<  (1. 

T  "  O 

■JLU 

329  C 

111,  J 

,  4  i  *T:iR  '  11 . 4’i  +AI1F  '  M,  J  ' 

330  C 

TO 

DFTEPI 1INE  SUPPORT  REAL 

TIONS 

33  1 

T  (  1 

'  =Ai1!1'  I  .  1 "'  -A  J 1 C  1 ) 

332 

T  ( 2 

)  =AI  111'  1 , 2  1  -  A  J  J  <  2) 

TT7 

IF'.; 

NMM.  EC! .  1  '  GO  TO  410 

34 


4  5  8 


339 

405 

340 

410 

341 

342 

343 

344 

415 

345 

C 

346 

C 

347 

348 

50O 

349 

350 

505 

354 

355 

356 


J  Ju 

359 


NtiriL  =1(1  in- 1 

NIK  IT =2+NMI1 

DO  405  11=  1 .  NMML 

DO  405  I=3,N11MT,2 

Tt  1 1  =h(  ir  i’.  n,  3  '  +411(141+1 .  i  >  —  h  j  j  <  1 1 


T (  2+tlMM+R =411(14(1111,  4  '  -AJJ  C2*NMH+2  > 
DO  415  1=1. HSS 
K  =  JNSS  < n 


PRINT  OUTPUT 
IJRITE  (6. 500) 

FORMAT'  ' 1' • 32X, ' MATRIX  ANALYSIS  OF  CONTINUOUS  BEAM',////) 

UP  I TE  > 6 -505i 

FORMAT ( IX. ' MEMBER ’ .5X, 'LENGTH' ,9X, *E' »9X. ' I' ,4X. 'LEFT  SHEAR' , 
1EFT  MOMENT ' .3X. 'RIGHT  SHEAR ',3X. 'RIGHT  MOMENT',//) 

DO  515  fl=  1  .  NMM 

UP  I  TE  f  6  •  5 1 0  i  II, XL  ( ( 1 J ,  E  <  M) ,  X I  <  11) ,  >'  AMI  K  M,  J ) ,  J  =  1 . 4) 

510  FORMAT (  L  :. I3.8X.F5.O.6X,F9.0. IX, FT. 1,4X,F7. 1.5X.F8. 1 .8X.F7. 1 . 

1.1.' 

515  CONTINUE 

UP  I TE ' 6 , 520 ) 

520  FORMAT' -. IX. ’DEGREE  OF  FREEDOM’ .21X, 'VERTICAL  DISPLACEMENT'  .! 
NOTATION- .  ■  i 


360 

IIP  I TE  >6.525) 

36  1 

tr  ~i  tr 
Jl.i 

FORMAT' 42X, "LEFT' 

362 

UR  I TE  >  6 . 530 ' 

363 

530 

FORMAT': 43X.  'OF-  .9 

364 

UP  I TE '6.535< 

365 

535 

FORMAT' 41X. ‘ JOINT 

366 

1  =  1 

367 

UP  I TE ' 6 . 540 ' I . DL ' 

368 

540 

FOPMhTi  7X.  13.38:-:. 

369 

M  l  J  =NM1 1-  1 

370 

DO  550  (1=1.11  JJ 

371 

1  =  1+2 

372 

UR  I TE ' 6 . 545 ' I . DP ' 

)  +SI1C 


3X,  ’  L 


6X,  F8 


IX.  'P 


13 


5b5 

573 


F OP t thi T i  7X,  13,26 X, F  1 Q . 5 . 2X.F10.5 ) 
LOUT [HUE 
I -1+2 

IJF'ITF >' (5,555>  I .  DR  ( HMM  -  3) 

F  iji- !  IhT i  77! -  13. 26X,  F  10.5) 

HP i fE ' 6.5601  I • PL f 1.2) 

FOFIIhT'7::.  I3.75X.F10.5) 

I'O  570  H=1.MJJ 
I  - 1  +2 
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Appendix  B 

Example  Output 


MATRIX  ANALYSIS  or  CONTINUOUS  MAM 
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A  computer  analysis  method  and  input  user  guideline 
to  determine  internal  forces,  reactions  and  deflections 
of  continuous  beams  with  rotational  slip  at  supports  is 
presented.  The  method  and  computer  program  are  applicable 
to  any  continuous  beam  structure,  although  developed 
specifically  for  two-span  floor  joist  analysis  and  design. 
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